Existence of global strong solutions for the shallow-water 
equations with large initial data 

Boris Haspot * 
Abstract 

This work is devoted to the study of a viscous shallow-water system with friction 
and capillarity term. We prove in this paper the existence of global strong solutions 
for this system with some choice of large initial data when > 2 in critical spaces for 
^ , the scaling of the equations. More precisely, we introduce as in JJJ a new unknown, a 

' effective velocity v ~ u + /iV In h {u is the classical velocity and h the depth variation 

of the fluid) with the viscosity coefficient which simplifies the system and allow 
us to cancel out the coupling between the velocity u and the depth variation h. We 
obtain then the existence of global strong solution if mo = hoVQ is small in 
' and (/iQ — 1) large in i?2^i- particular it implies that the classical momentum 
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itIq — HqUq can be large in , but small when we project uiq on the divergence 

. field. These solutions are in some sense purely compressible. We would like to point 

out that the friction term term has a fundamental role in our work inasmuch as 
coupling with the pressure term it creates a damping effect on the effective velocity. 
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1 Introduction 



We consider the viscous shallow water model with friction and capillarity term. This 
model is also called by the french community the Saint- Venant equations and is generally 
I used in oceanography. Indeed it allows to model vertically averaged flows in terms of the 

' horizontal mean velocity field u and the depth variation h due to the free surface. In the 

rotating framework, the model is described by the following system: 



d 

—h + div(/in) = 0, 

d Vh ^ ' 

— (hu) + div( hu u) — div(/ihDu) + —r + rhu = divK, 
ot Fr^ 

where divK is the free surface tension tensor which reads as follows: 

divK = V{hK{h)Ah + + Hk {h))\Vh\^) - dw{K{h)Vh (g) Vh) . (1.2) 



a 



K is the coefficient of free surface tension and is a regular function of the form K{h) = Kh 
with a G R. In the sequel we will assume that a = —1. Fr > denotes the Froude 
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number. System (jl.ip is supplemented with initial conditions 

h/t=Q = ho, (M/t=o = "^0- (1-3) 

This model is derived from the three-dimensional Navier-Stokes equations with free sur- 
face, where the normal stress is determined from the air pressure and capillary effects. 
The turbulent regime (r > 0) is obtained from the friction condition on the bottom, see 
[23]. n is the viscosity coefficient and verifies ^ > and Du = (Vn+* Vn)/2 is the strain 
tensor. 

Several physical models arise as a particular case of system (jl.ip : 

• when K = r = 0, (jl.ip represents compressible Navier-Stokes model with shallow- 
water viscosity coefficients. 

• when K > and r = 0, then (jl.ip describes the Korteweg system which models 
mixture liquid-vapor. 

We would like to point out also the theoretical aspect of system (jl.ip . Indeed in the 
case r = the system (jl.ip corresponds to the classical Korteweg system which models a 
liquid-vapour mixture. Let us mention that the Korteweg system also is used in a purely 
theoretical interest consisting in the selection of the physically relevant solutions of the 
Euler model by a vanishing capillarity-viscosity limit (in particular when the system is 
not strictly hyperbolic which is typically the case when the pressure is Van der Waals). 
Indeed in this case at least when A'^ = 1 it is not possible to apply the classical theory 
of Lax for the Riemann problem (see [20]) and of Glimm (see [13] ) with small initial 
data in BV in order to obtain the existence of global entropic solution (we refer also 
to the work of Bianchini and Bressan see [2] for the uniqueness). In particular in this 
direction, recently in [5] with F. Charve we prove that the global strong solution of the 
Korteweg system in one dimension (we obtain also in this paper the existence of global 
strong solution in one dimension for Korteweg system) converges in the setting of a 7 
law for the pressure {P{p) = ap^ , 7 > 1) to entropic solution of the compressible Euler 
equations. In particular it justifies that the Korteweg system is suitable for selecting the 
physical solutions in the case where the Euler system is strictly hyperbolic. The problem 
remains however open for a Van der Waals pressure. 

Now before investigating the problem of global strong solution for the system (jl.ip . we 
would like to recall the energy inequalities associated to this system and in particular 
describing some results about the existence of global weak solutions for the system (jl.ip . 
Let /i > be a constant depth variation (in the sequel we will assume that h = 1), and 
let n be defined by: 

so that -P(s) = sll (s) — n(s) , n [h) = 0. Multiplying the equation of momentum 
conservation in the system (jl.ip by u and integrating by parts over R^, we obtain the 
following estimate: 

/ [\h\u\'' + (n(/i) - n(/i)) + ^|v/i|2) {t)dx + \ f I ph\D{u)\^dxdt 

- L (^ + (^(M-^(M) + ^|v/^o^dx, 
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Here we can observe that if we assume that the initial data are such that: 



ho H—) e L\ y^V/io €L , 

e (1.5) 

y%\uo\ G L\R^), 
then we have the following estimates: 

h\u\^ G L°°(Li), (n(/i) - U{h)) e L°°(L^), 

^|V/i|2 G L~(Li), Z^^x G L\L2). ^^'^^ 

One of the main difficulty in order to obtain the existence of global weak solution consists 
in dealing with the quadratic terms of the capillarity tensor. Indeed in order to have the 
stability of a sequence {hn,Un)neN of global weak solution for the system (jl.ip . it is 
crucial to give a sense to the quadratic terms in the gradient of the depth variation 
K{hn)'Vhn <8) V/in which are only uniformly bounded in L°°(L^(]R^)). In particular it 
implies only a convergence in the sense of the measure of K{hn)'Vhn ^ V/i„ which is not 
sufficient to conclude by standard compactness argument. That is why the problem of 
the existence of global weak solution remains open in the general case. Before stating 
our main result on the existence of global strong solution, we would like to recall what 
is known on the existence of global weak and strong solution in the different standard 
configurations. And in particular we would like to emphasize on the results obtained in 
|14^ [T5] where we have discovered new entropies on the Korteweg system (r = 0) for a 
specific choice on the capillarity {k{p) = ^ with k > 0), allowing in particular to prove 
the existence of global weak solution. It is precisely with this type of capillarity that we 
will work with in the sequel. 
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Existence of global weak solutions 

Case K = and r = 0, the compressible Navier-Stokes system 

When the viscosity coefficients are constant and the pressure is a 7 law P{h) = ah'^ , with 
a > and 7 > 1, P-L. Lions in [21] proved the global existence of variational solutions 
{h, u) to (mi with K = r = 0for7>fif7V>4, 7> ^ if N = 2,3 and initial data 
{po,mo) such that: 

I 1 2 

u{ho)-uCh), ^glHm"^). 

tio 

These solutions are weak solutions in the classical sense for the equation of mass conser- 
vation and for the equation of the momentum. Notice that the main difficulty for proving 
Lions' theorem consists in exhibiting strong compactness properties of the height h in 
L^g^ spaces required to pass to the limit in the pressure term P{h) = ah'^. 
Let us mention that Feireisl in [10] generalized the result to 7 > ^ i'^ establishing that we 
can obtain renormalized solution without imposing that h G L^g^., for this he introduces 
the concept of oscillation defect measure evaluating the lost of compactness. 
Concerning the shallow- water system when = /i/i, the main difficulty when dealing 
with vanishing viscosity coefficients on vacuum is that the velocity cannot even be defined 
when the density vanishes. In particular we lose the information on Vu in L'^{{0, T) xR^). 
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The main difficulty, to prove the stabihty of the solutions, is to pass to the limit in the 
term pu®u (which requires the strong convergence of ^/pu) . Mellet and Vasseur in [22j 
obtain the stability of global weak solution by using new entropies on the velocity and 
the density. 

Case r = 0, Korteweg system 

In the capillary case in contrast to the non capillary case, we can easily deal with the 
pressure term in order to obtain stability results. However let us emphasize at this point 
that the energy estimates do not provide any L°° control on the density from below or 
from above. Indeed, even in dimension N = 2, functions are not necessarily locally 
bounded. Thus, vacuum patches are likely to form in the fluid in spite of the presence of 
capillary forces, which are expected to smooth out the density. It explains why it is so 
difficult to obtaining the existence of global strong solution in dimension N = 2 and even 
global weak solution. In [IT], we obtain the existence of global weak solution for specific 
choices of the capillary coefficients and with general viscosity coefficient but with small 
initial data in the energy space. More recently in |14] and |15] . we prove the existence 
of global weak solution with large initial data for the Korteweg system (when K{h) = j^) 
with Saint- Venant viscosity coefficients. Indeed with this choice of capillarity coefficient 
we are able to obtain new entropy inequalities. To do this, we introduce a new unknown 
V = u+ /iV In p called effective velocity, and we are able to show some gain of integrability 
on V which allow us to deal with the terms of the type pu® v (we refer to [15J for more 
details). 

Existence of global strong solutions 

Case K = and r = 0, the compressible Navier-Stokes system 

We refer to |18] for the existence of global strong solution with small initial data in critical 

JV ——I 

space for the scaling of the equations. More precisely (/iq — 1, uq) belongs to B^^ x B^^ ^ 
for suitable chosen p and pi . 

Case r = 0, Korteweg system 

Let us mention briefly that the existence of global strong solutions in critical spaces for 
the scaling of the equations for > 2 is known since the works by R. Danchin and B. 
Desjardins (when the viscosity coefficient and the capillary coefficients are constant, see 

JV JY._i 

also |16] ) in [8] where the initial data (/iq — 1, /lo^o) belong to the Besov spaces B21 '>^B^^ 
and are chosen small enough. In [14j, we improve this result by working in a larger space 

JV i^ — 1 

of initial data, more precisely (/iq — 1, /io^o) is in the Besov space -62^00 ^ -^2^00 • ^^^^ 
paper for the same reason than in [15], we are working with specific capillary coefficient 
K{h) = because we are able to exhibit a specific structure on the unknown In h. In 
particular we are able to work with vortex initial data on the initial velocity in dimension 
N = 2. 
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1.1 Derivation of the model and non trivial explicit solutions 

The choice of system (jl.ip is motivated by its energetic consistency, which has been 
stressed out from a physical point of view in [11] . Compared with the Korteweg system ( 
see [19]) we also take into account the friction term rhu physically justified to model the 
friction condition on the bottom of the ocean ( see [12] )■ Naturally on a mathematical 
point of view, the term rhu does not add any difficulties for obtaining global weak 
solution or global strong solution with small initial data, however in our study this term 
shall turn out to be essential to ensure the existence of global strong solution with large 
initial data. Indeed coupled with the pressure term, he shall introduce a damping effect 
on a new unknown called effective velocity and introduced in |14t [TC] . Roughly speaking 
this friction term allows to cancel out the coupling between the height h and this effective 
velocity it is one of the main difficulty in order to obtain the existence of global strong 
solution for compressible Navier-Stokes equation (see [HEKTS], in this case the coupling 
is between the velocity and the density). Indeed in this last case, it is difficult to obtain 
a damping effect on the density and the coupling between the velocity and the pressure 
terms impose a smallness condition on the initial density. In our case the fact that with 
the friction term we can rewrite the pressure term in some sense as a velocity, more 
precisely as the effective velocity allows to avoid any condition of smallness on the initial 
depth variation /iq. We shall come back more in details on these considerations in the 
proof of theorem 11.11 we now would like to introduce this notion of effective velocity 
which is one crucial tool for the proof of our result. 

More precisely we are going to explain why with specific choice on the capillarity, Froude 
and friction coefficient (as in |14] and [15j) we can exhibit a new structure onthe system 
via the introduction of this effective velocity v. In particular in |14^ [T5] we obtain new 
entropies which allows us toprove the existence of global weak solutions. We are now 
considering in the sequel the following physical coefficients: 

K 1 
K.{h) = -, K = if^ and — = r/i, (1.7) 
h rr'^ 

with > 0. By computation (see PH])) we obtain the simplified system: 

f dth + d\Y{hv) - fiAh = 0, 

1 hdtv + hu • Vv — div(^h Vv) + r h v = 0, 

with V = u + (xV In h the effective velocity. For more details on the computation, we 
refer to [15]. When we write the system (jl.Sp in function of the momentum m = hv, the 
system reads as follows: 

{dfh + divm — fiAh = 0, 
m (1.9) 
dtm + div(— (g) m) — //Am + r m = 0, 
h 

In particular we can observe that m = and hi such that: 

{dthi — riAhi = 0, 
hi{0,x) = hoix). 
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is a particular solution of (jl.lOp . When we consider the system (jl.ip . it means that 
(/ii,Mi = /xVln/ii) is a non trivial solution. In particular this solution is purely com- 
pressible as curlui = 0. 

In the sequel we are going to work around this non trivial solution (/ii, 0) for the system 
p.lOp . If we now consider the system (jl.lOp . we can observe that the coupling between 
the momentum m and the height h via the pressure term has disappeared. As we ex- 
plained previously the friction term and the pressure introduce here a damping effect via 
the new term rm. It is now possible to obtain the existence of global strong solution 
when we assume only a condition of smallness on mo in order to deal with the non linear 
term div(^ (8) m) and in particular to use smallness argument for the term ^. 



Results 

Our main motivation concerns the existence of global strong solution with large initial 
data on the height ho but also on the initial velocity uq. More precisely we will obtain 
global strong solution with a family of initial velocity uq such that the projection on 

gradient vector field has large norm in ■ 

We prove global well-posedness for system ()1.10l ) in critical Besov space. To do this, 
we shall work around a constant depth variation h = 1, and to do this we introduce the 
following definition. 

Definition 1.1 We set: 

qo = ho - 1. 

We can then rewrite the system (|1.10p in function of {q,m), it gives: 

{dtq + divm — fiAq = 0, 
5tm + div(— (8) m) — /iAm + r m = 0, ^ ^ 

h 

Theorem 1.1 Suppose that we are under the conditions |./.7y . Assume that mo G 

JV 

and qo G with Hq > c > 0. Then there exists a constant eq depending on ^ such 
that if: 

||"io|| 4-1 < ^0, 

^2,1 

then there exists a unique global solution {q, m) for system with h bounded away 
from zero and, 

~ K —A-1 ~ — — 1 ——1 — -1-1 

h £ C(M+, B2\) n L^(M+, 1 ) and m G C(M+; ^s'l ) n L^(M+, n B^^ ). 



Remark 1 In this theorem, for the first time up my knowledge we obtain a result on the 
existence of global strong solution for a compressible system without assuming any small- 
ness hypothesis on the density as it is classically the case (see jQ, [71 [3 ES'J- Furthermore 

the physical moment uiq = HqUq is large in B21 , indeed we have: 

m'Q = mo- AiVpo, 
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^-1 

and here ruQ is small in B^i hut Vpo could be arbitrary large. It means that for this 
family of large initial data ( when vtlq is the sum of a small momentum and of the gradient 
of the density), we have the existence of global strong solution. In particular the projection 

on the divergence field issmall in but the projection on the gradient field is large 

in . It is essentially related with the compressibility of the system. 

Remark 2 In the same way than in fl^}, we could consider the unknown (In /i, u) and 

in this case we could obtain the existence of global strong solution with small initial data 

in (^2^00' -^2^00 ) which is a larger space than in theorem In particular it allows 

us to deal with the problem of vortex initial data, indeed in dimension N = 2 we could 
choose initial velocity such that curluo is a bounded measure and in particular a Dirac 
Sq. However it would be not clear how to get global solution without smallness hypothesis 
on the initial density, indeed in the term u ■ Vv that we write under the form v ■ Vv — 
liVlnp ■ it is not clear how to deal with the quadratic term Vlnp • Vf . 

Remark 3 Our method may be adapted to the framework (that is we now consider 

JV — — 1 

go S Bp T^, niQ € Bp^^ ) when I < p < +00 as in f^. 

Our paper is structured as follows. In section [21 we give a few notation and briefly 
introduce the basic Fourier analysis techniques needed to prove our result. In section [3l 
we prove the theorem dl.il 



2 Littlewood-Paley theory and Besov spaces 

Throughout the paper, C stands for a constant whose exact meaning depends on the 
context. The notation A < B means that A < CB. For all Banach space X, we 
denote by C([0,T],X) the set of continuous functions on [0,T] with values in X. For 
p G [1, +c«], the notation L^'(0, T, X) or Vrp{X) stands for the set of measurable functions 
on (0,T) with values in X such that t — )• belongs to LP{0,T). Littlewood-Paley 

decomposition corresponds to a dyadic decomposition of the space in Fourier variables. 
We can use for instance any cp G C°°(M^), supported in C = G < ICI < §} such 

that: 

^v=(2-'e) = i if e/o. 

Denoting h = J^~^ip, we then define the dyadic blocks by: 



kU. 

k<l-l 



Aiu = ip{2-^D)u = 2^^ / h{2^y)u{x - y)dy and Siu = V A 
Formally, one can write that: 

u = ^ AfcU. 

This decomposition is called homogeneous Littlewood-Paley decomposition. Let us ob- 
serve that the above formal equality does not hold in S (R^) for two reasons: 
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1. The right hand-side does not necessarily converge in S 



s>N\ 



2. Even if it does, the equality is not always true in S (M^) (consider the case of the 
polynomials). 

2.1 Homogeneous Besov spaces and first properties 

Definition 2.2 For s G M, p G [1, +oo], q G [1, +oo], and u G 5'(M^) we set: 

/gz 

The Besov space B^^ is the set of temperate distribution u such that \\u\\b^^ < +oo. 

Remark 4 The above definition is a natural generalization of the nonhomogeneous Soholev 
and Holder spaces: one can show that B^ ^ is the nonhomogeneous Holder space and 
that B2 2 is the nonhomogeneous space . 



Proposition 2.1 The following properties holds: 

1. there exists a constant universal C such that: 
C-^\\u\\bs^^ < ||Vn||^.-i < C\\u\\bs^^. 

2. If Pi < P2 and ri < ra then B^^^^.^ ^ Bllr}'^'^^~^'^^\ 

' I I 

3. Bp j.^ ^ if s > s or if s = s and n < r. 

Let now recall a few product laws in Besov spaces coming directly from the paradifferen- 
tial calculus of J-M. Bony (see [3j) and rewrite on a generalized form in [?] by H. Abidi 
and M. Paicu (in this article the results are written in the case of homogeneous sapces 
but it can easily generalize for the nonhomogeneous Besov spaces). 

Proposition 2.2 We have the following laws of product: 

• For all s G M, {p,r) G [l,+oo]^ we have: 

\Wv\\b=p^, < C{\\u\\loo\\v\\b=^,. + \\v\\l^\\u\\b=^^^) . (2.12) 

• Let {p,Pi,P2,r,XiA2) G [l,+oo]2 such that:^ < ^ + ^, Pi < A2, P2 < Ai, i < 



^ + and - < 7^ + 1^. We have then the following inequalities: 
ifsi + S2 + Nm/{0, 1 - J- - ^) > 0, + f < f and S2 + < f then: 



\\uv\\ , , wr 1 -1- 1 i^5JI'u|Ir^i Iklln^a , (2.13) 

when si + g = ^ (resp s2 + ^ = ^)we replace \H\b;iJMb;i^ f^esp \\v\\b;i^) 
by \W\\b;iJHb;1, (resp \HB;i^nL^)' ^1 + = ^ (^^id S2 + ^ = ^ we take 
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r = 1. 

// + S2 = 0, si G (f - ^, ^ - f 1 and^ + ^ <1 then: 



_jv(i+i_i) < lltill^n 1^11^=2 . (2.14) 



^Pl P2 P Pl,l P2' 

If \s\ < ^ for p > 2 and —^■ < s < ^ else, we have: 



\uv\\b'> < C'IIuIIrs N . (2.15) 

I II o.r — II II o.rii II "TT ^ ' 



Remark 5 In the sequel p will he either pi or p2 and in this case j = — if pi <P2, 



'^esp i = ^ - ^ ifP2<Pi. 
Corollary 1 Let r £ [1, +00], 1 < p < pi < +00 and s such that: 

N_ 

then we have if u £ Bp ,, and v G -Bpi\oo H L°° : 

WuvWb" < ClliillB" \\v\\ JV 



A pi P2 



p,r — II II p,r ■■ ■■ 



RPi nr,°° 

-^pj^jOol 1-^ 



The study of non stationary PDE's requires space of type L''(0,T, X) for appropriate 
Banach spaces X. In our case, we expect X to be a Besov space, so that it is natural to 
locahze the equation through Littlewood-Payley decomposition. But, in doing so, we ob- 
tain bounds in spaces which are not type LP{0, T, X) (except if r = p). We are now going 
to define the spaces of Chemin-Lerner in which we will work, which are a refinement of the 
spaces Lj.{Bp .i.). 

Definition 2.3 Let p G [l,+oo], T G [l,+oo] and si G M. We set: 

We then define the space L^{Bp^^) as the set of temperate distribution u over (0, T) x 
such that \\u\\jp (ns-i^ , < -|-cxd. 

We set CriBp^r) = -^r(-^p,V) ^ C([0, T], i?p\,). Let us emphasize that, according to 
Minkowski inequality, we have: 



Remark 6 It is easy to generalize proposition \2.S\. to Lj,{Bp^^) spaces. The indices si, 
p, r behave just as in the stationary case whereas the time exponent p behaves according 
to Holder inequality. 
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In the sequel we will need of composition lemma in L^(i?p spaces. 
Lemma 1 Let s > 0, (p,r) G [1, +oo] and u G L^{B^^r) ^ L^{L°°)- 

1. Let F G W;W+^'°°(R^) such that F(0) = 0. Then F{u) £ Lf^{B^^^). More precisely 
there exists a function C depending only on s, p, r, N and F such that: 

2. Ifv,u€ L^iB^^r) n L^{L°°) and G G w\'}^'^'°° {W^) then G{u) - G{v) belongs to 
L?p{Bp) and there exists a constant C depending only ofs,p,N andG such that: 

\\Giu) - ^(t;)!!^^ < C(||n||ioo(^^), \\v\\l^^l^)){\\v - u\\^^^^^ Jl + \\u\\l^^l^) 

+ II^IIl-(l-)) + 11^^ - + 

Now we give some result on the behavior of the Besov spaces via some pseudodifferential 
operator (see [1]). 

Definition 2.4 Let m G M. A smooth function function f : — )■ M zs said to be a 5™ 
multiplier if for all muti-index a, there exists a constant Ga such that: 

vcgm^, <c„(i + ieir-i"i. 



Proposition 2.3 Let m G M and f he a 5™ multiplier. Then for all s €z M and 1 < 
p,r < +00 the operator f{D) is continuous from Bp ^. to Bp~"^. 

Let us now give some estimates for the heat equation: 

Proposition 2.4 Let s G M, (p, G [li+oo]^ and 1 < /?2 < Pi < +oo. Assume that 
uq G Bpj. and f G L^ {Bp^r^^"^^ '^^) . Let u be a solution of: 

dtu — fj.Au = f 
Ut=0 = uo . 

Then there exists C > depending only on N, /i, pi and p2 such that: 
If in addition r is finite then u belongs to C([0,T], Bp ^.). 
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3 Proof of the theorem 11.11 

The existence part of the theorem is proved by an iterative method. We define a sequence 
(g",m") such that: 

dt(f — fiAq^ + rmdivm^ = 0, 
dtrmP — fiAm^ + rm^ = 0, 
(g°,m°) = {qo,mo). 
Assuming that [q^^m^) is in Et with: 

Et = {Ct{b,'^,) n lUb^/ )) X {CT{Bi, ) n L^Bi/ n B^ - ))^, 

we define then q^^^ = q^ + Q^'^^j mP"^^ = + m^^^ such that (g„+i,m„+i) be the 
solution of the following system: 

dtq^+^ - nAq^+^ + divm"+^ = 0, 

[(r+\m"+i)/i=o = (0,0), 
with: 

= - div(— m'^) 

We also set: /i" = + 1. 

1) First step, uniform bounds: 

Let e be a small parameter and by proposition 12. 4| we have for any T > 0: 
llo'^ll jv jv , 9 < C||(7oll N , 

Tl 2.1> TV 2,1 2,1 

||m II jv_-^ jv_i iv_|_i < C||mo|| N_-^. 

^t(^2,1 )^-^t(-'^2,1 ^-^^2,1 ) -^2,1 

We are going to show by induction that for e > small enough: 
(Vn) ||r,m")||^, <e. 

As {qo,fhQ) = (0,0) the result is true for n = 0. We suppose now (Vn) true and we are 
going to show (Vn+i)- 

To begin with we are going to show that 1 + is positive. Indeed we have: h? = h\-\- 
such that: 

dth\ - nAhl = 0, 
{hi)/t=o = ho. 

and: 

dthl - fiAhl = -divm°, 
(/ii)/t=o = ho- 

By proposition ()2.4p and ()3.16p we have for any T > 0: 

||/i^|| N < C\\mo\\ N (3.17) 
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By maximum principle, we have for any t > 0: 

h^(t,x) > min hnix) > c > 0. 

We deduce that for tj = ||mo|| n (at least inferior to ^ with the C of (|3.17p ) small 
enough and any t > 0: 

hyt,x)/geqj >0, 

and 

g°(t,x)>f -1. 
and by definition of and the assumption nP„ that: 

9"(t,x)>--l-e. 

In particular for e small enough at least e < |, we deduce that: 

/i*^ = 1 + > £ > 0. (3.18) 

In order to bound (g", m") in Et, we shall use proposition l'i . 41 and in particular estimating 
in L}p{B2i ). By using proposition 12.21 (jS.lSp and lemma[Tl we obtain: 

||div(— (g)m'^)|| N , <||— ®m''|| n 

LliB^,) LliB^y L^iBZ) L^iB^i) 

Therefore by using p.l9p . the proposition 12.41 and (Vn) we obtain for any T > 0: 

||(r+\m"+i)||F. <C(||mOf ^ +ef{e+\\q'\\ . 

By choosing rj = e and e < 2c(2+||g"|| ^^^^ implies {V)n+i- We have shown by 

induction that {q^,m^) is uniformly bounded in Ft for any T > 0. 



(3.19) 



Second Step: Convergence of the sequence 

We shall prove that {q^,m^) is a Cauchy sequence in the Banach space Ft, hence con- 
verges to some {q,in) E Ft- 
Let: 

= - and 5m" = m"+i - m". 
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The system verified by {6q^,5m"') reads: 

' dt5q'^ - ^A(5g" + div(5m" = 0, 
< atJm" - ^A(5m" + r5m" = - Gn - 1, 
_ 5g"(0) = , (5n"(0) = 0. 

Applying propositions 12.41 and using (Vn), we get for any T > 0: 

By using proposition 12.21 and lemma [U we get: 

\\{5q^,5m^)\\F^ < Ce\\{6q'''\6m^-')\\F^. 

So by taking e enough small we have proved that {q^,ra^) is a Cauchy sequence in Ft 
which is a Banach space. It implies that {q'^,m^) converge to a limit {q,m) in Ft- It is 
easy to verify that {q,m) is a solution of the system (jl.lip . 



3) Uniqueness of the solution: 

The proof is similar to the proof of contraction, indeed we need the same type of estimates. 
Let us consider two solutions in Ft- {qi,mi) and {q2,m2) of the system (jl.lip with the 
same initial data. With no loss of generality, one can assume that (gi, mi) is the solution 
found in the previous section. We thus have: 

(n) qi{t,x)>-^. 
We note: 

5q = q2- qi, 6m = m2 - mi, 

which verifies the system: 

dt5q — fiASq + divJm = 0, 

dtdm — fiA6m + r6m = — div(— ^ (8) mi) + div(— ^ mi) 

hi hi 

By using proposition 12.21 12.41 and lemma [1] on [0, Ti] with < T we have: 

\\{dq,6m)\\E^ < A{T)\\{6q,6m)\\E^, 

such that for T small enough A{T) < ^. We thus obtain: 6q = 0, 6m = on [0, T]. And 
we repeat the argument in order to prove that: 6q = 0, 6m = on This conclude 
the proof of theorem 11.11 □ 
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